The integrated density of states of the random graph Laplacian 
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We analyse the density of states of the random graph Laplacian in the percolating regime. A 
symmetry argument and knowledge of the density of states in the nonpercolating regime allows us 
to isolate the density of states of the percolating cluster (DSPC) alone, thereby eliminating trivially 
localised states due to finite subgraphs. We derive a nonlinear integral equation for the integrated 
DSPC and solve it with a population dynamics algorithm. We discuss the possible existence of a 
mobility edge and give strong evidence for the existence of discrete eigenvalues in the whole range 
of the spectrum. 



The eigenvalue spectrum of sparse random matri- 
ces is a fascinating and largely unsolved problem with 
widespread applications ranging from transport in disor- 
dered systems, graph theory and optimization problems 
to nuclear physics and QCD [TJ[2]. In this paper we con- 
sider a prototype of such a random matrix, the Laplace 
operator on a mean-field random graph with N nodes, i.e. 
a graph where a link between two arbitrary sites is either 
present with probability p = 2c/ N or not present with 
probability 1 — p. The constant 2c is the mean connec- 
tivity of the graph. The Laplace operator on this graph 
is a matrix where for i j Tjj = — 1 if the nodes i 

otherwise, 



with n nodes is given in the macroscopic limit by [lOj 



Y^jjti^ij- The entries on 



and j of the graph are connected and T^- 
while on the diagonal Yu - 

the diagonal are thus correlated to the random entries 
outside the diagonal. 

Even though the computation of the density of states 
for a mean-field random graph has been reduced to an in- 
tegral equation [3HS] , a complete solution is still missing. 
In the limit of infinite coordination, c — > oo, Wigner's 
semi-circle law is recovered. For any finite c, Lifshitz tails 
were shown to exist in the integrated density of states [6] . 
Beyond these asymptotic results a number of approxima- 
tions have been used to compute the spectrum approx- 
imately, such as effective medium theory, single defect 
approximation, moment expansions and numerical diag- 
onalisation [THS]. 

In this paper we show first that the density of states 
of the percolating cluster (DSPC) can be isolated using 
a symmetry argument and our knowledge of the density 
of states below the percolation threshold. Second, we 
derive an integral equation for the integrated density of 
states which can be solved reliably with a population 
algorithm. The numerical solution reveals jumps in the 
integrated DSPC in the whole range of the spectrum, 
calling in question the existence of a mobility edge. 

Model and symmetry: The model shows a percola- 
tion transition at c cr u — 1/2. Below this concentration 
there is no macroscopic cluster and almost all finite clus- 
ters are trees. The average number of tree clusters T n 
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In particular the total number of clusters per particle is 
T to t(2c) = 1 — c. For c < 1/2, the spectrum consists of 
a very complicated, but countable set of <5-peaks which 
can be calculated iteratively [IT]. Above the percola- 
tion threshold c > 1/2 a percolating cluster coexists with 
many finite clusters, which are also trees. The fraction 
of sites in the macroscopic cluster, Q(c), is the solution 
of 1 — Q(c) = exp(— 2cQ(c)). Alternatively we rewrite 
Q(c) = 1 — and obtain x(c) as the solution of 



x(c) 



-x(c) 



2ce 



-2c 



(2) 



This equation has two solutions, a trivial one with x(c) = 
2c and a nontrivial one with x(c) — 2c* such that c* > \ 
if c < \ and vica versa. This nontrivial solution allows 
one to establish a symmetry for the number of trees above 
and below the percolation threshold. Using Eq. ([2]), we 
can rewrite Eq. ([I]) according to 

r n (x(c)) = 4V"( 2c ) or T «( 2c *) = 4 T »( 2c )- ( 3 ) 
x(c) c* 

Hence the number of trees above the percolation thresh- 
old is simply related to the number of trees below the 
percolation threshold. 

Density of states This relation allows us to com- 
pute the density of states of the percolating cluster alone, 
which is the quantity of primary interest. It is known 
that the density of states, -D(fi), of the infinite cluster 
contains at least some S peaks, so that a population 
dynamics algorithm [5] cannot be applied. In this pa- 
per we instead compute the integrated density of states, 
A(Q) = dQ,'D(Q,'), which according to measure the- 
ory |12) . may be decomposed into a singular part and 
an absolutely continuous part. The absolutely contiuous 
part may itself consist of two contributions, eigenvalues 
stemming from localised eigenvectors which happen to lie 



2 



dense and each have vanishing weight in the thermody- 
namic limit, and the eigenvalues stemming from nonlo- 
calised eigenvectors. We may thus write 

A(Q) = Ai OC;disc (f7) + A loc , cont (ft) + A nonloc (fi). (4) 

Often in random matrix problems there is a mobility 
edge, i.e. a value fio such that all eigenvectors corre- 
sponding to eigenvalues f2 < f2 are localised and all 
eigenvectors corresponding to Q > fio are nonlocalised 
(or vice versa) . We will show here that such a sharp edge 
does not seem to exist in our problem as we find dis- 
crete eigenvalues even in the region where nonlocalised 
eigenvectors lie. This shows that a naive approach using 
the inverse participation ratio in order to find the on- 
set of nonlocalised eigenvectors will not work since the 
localised eigenvectors do not disappear where the nonlo- 
calised ones start, so the inverse participation ration will 
not drop down to as it would at a true mobility edge. 

The density of eigenvalues, {Oj}^, of the Laplacian 
matrix T is defined by 



1 N 1 

D(n,c)= hm — ys(n-fii)= iim — TW(fi-r) 

1=1 

(5) 

Here • denotes the average over all realizations of connec- 
tivity for a given c. To compute the density of eigenvalues 
we introduce the resolvent 



G(fi, c) 



lim — Tr - 



jv-x» n r - n 



(6) 



for complex argument Q — 7 + ie, e > 0. In the limit 
e — > 0, we recover the spectrum from the imaginary part 
of the resolvent according to 



D(Q,c) = ~ lim SG(fi + ie,c). 

7T 40 



(7) 



In the percolating regime c > I the macroscopic cluster 
coexists with many finite ones. In order to study localised 
states of the macroscopic cluster it is essential to isolate 
the density of states of the macroscopic cluster only. We 
use Eqs. ^ and ^ to decompose the resolvent into two 
contributions, one from the percolating cluster and one 
from the finite clusters: 



G(Sl, c) = G pcTC (n, c) + —G(Q, c* 
c 



(8) 



The above relation allows us to compute the density of 
states of the percolating cluster alone from the full den- 
sity of states above, G(tt,c), and the density of states 
below the percolation threshold, G(il,c*) for c* < i, 
which is known 1111. 



The average over all realizations of connectivity is per- 
formed with the replica trick, resulting in a nonlinear 
integral equation for g c ,n{p) (cf. Eqs. (16) and (17) in 
Ref. 0) 



9c,n{p) 
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with <? c ,fi(0) = 2c. Here l v (z) are the modified Bessel 
functions of the first kind. The solution of Eq. ([9| yields 
the resolvent [3] 



G(fi,c) = -1 + - / 
2c J 



dppg c ,n(p) ( 10 ) 



Nonpercolating regime: The analytical solution of 
Eq. j9) for c < 1/2 was given in [llj . Briefly, it is 

/oo ^ 
da Cj n(A)e"5 Ap 

(11) 

with coefficients and Zk to be defined below. The sum 
can also be expressed as a Riemann-Stieltjes integral with 
weight function a c ,u(A) = J2k a k&(^ — z k) {&{■ ■ ■) is the 
Heaviside function). This formulation will be useful later. 



The coefficients and can be grouped in infinitely 
many "classes." The coefficients of class n + 1 are given 
recursively by the relations 



n (2c4 n) )'» 
11 u 



(n) 



n) 



(12) 

(13) 



Class contains 
u 



e 2c and z^ } - r — 



(the upper index denotes the class) 
only one element, namely a Q ^ 
The index on the left hand side is a finite sequence (Ik) 
of nonnegative integers. In order to proceed to the next 
stage of the recursion, (Ik) must be mapped to a natural 
number m since for the calculation of, say, 



in) 

necessary to access the coefficients z m of the previous 



it is 
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iteration. Such a mapping is possible because the set of 
sequences {{Ik)} is countably infinite. It was shown in 
[TT] that the correct way to do the mapping is to choose 
m = lkM k and to let M (formally) tend to infinity at 
an appropriate point. Note that class n + 1 also contains 
all coefficients from class n. See QT] for details. 

These classes give us an infinite hierarchy of coeffi- 
cients, each recursion step adding infinitely many coef- 
ficients to the previous ones. Note that the coefficients 
in class n constructed in this way are not an approxima- 
tion but constitute (part of) the exact solution of Eq. ^ . 
Only the total weight of coefficients J^k °fc f a hs short of 
1 when stopping the recursion at a finite n. This solu- 
tion of Eq. ([9| leads to a density of states consisting of (5 
peaks which are located at those Cl where = 0. 

Percolating regime: In complete analogy to the re- 



solvent, we can decompose g c ,n(p) = 9^n(p) + 9c,n(p) 
into a part <7?q C (/£>) pertaining to the infinite cluster and 
a part which is equal to the (known) solution g c *.n(p) 
of the integral equation at c*. This decomposition has 
the advantage that we can directly obtain the density of 
states of the infinite cluster by deriving an equation for 



9c o C (/°)' which m analogy to Eq. (11 1 can be represented 
as 



g%?(p) = (2c -2c*) 



d& Cl n(A) e 



(14) 



The function fc c .n(A) is normalized such that 
J^°oo ^c,n(A) = 1 and the "initial condition" g c ,n(0) = 2c 
is being taken care of by the prefactor 2c — 2c* . This 
ansatz is plugged into the integral Eq. ^ to yield 



b c>n (A) = 2c*^a„ £ 



n=0 M=\ 



(2c -2c*) 

mT 



M-i r°o 



<#>c,n(Ai) 



<#> c ,o(Am) x 8 A- 



1 - 



(15) 



This equation can be solved numerically by running a 
population dynamics algorithm for the coefficients Zk at 
c* below the critical point in parallel to a population 
dynamics for a A-population at c above the critical point, 



for which the Zk and their weights arc needed as input. 

Given a c * t a(X) and 6 c ,n(A) the integrated DSPC can 
be computed according to (see [T3"]): 



2A^ c (Cl,c) = (l + J°° db Ct n{\) ^gB( I ^ I )+c Sgn (A-l)) -(c-c*) j° 



d6 c ,n(A)d6c,fi(A')sgn( 



A-l 



r 



A') 



< < I da c ., n (A)sgn( x ^ T )-c* / (da c , >u (X)db c>u (X') + db c , Q (\)da c *, Q (\'))sgn( ^ - A') j . 



Discussion: We have developed a systematic ap- 
proach to compute the integrated DSPC. We stress that 
the (nonintegrated) DSPC could not be reliably obtained 
from a population dynamics algorithm. Naively it is 
given in terms of & c ,fi(A) by 



D pcrc (n,c) 



1 



dpp d& C; o(A)e 2 



db c . n (X)S(X) 



-IV 



(17) 



The above density of states would be very simple if 
b c ,n(X) was differentiable with respect to A, as it would 
then be equal to (fe c ,n)'(0). It is however known that 
the density of states of the percolating cluster contains 
5 peaks also for c > 1/2. Hence the position and weight 
of the S peaks can not reliably be obtained from a popu- 



lation dynamics algorithm - a problem already encoun- 
tered in the nonpercolating regime. It is thus essential to 
obtain the integrated density of states directly from pop- 
ulation dynamics without going to the density of states 



first and integrating, and this is what Eqs. (151 and (16) 
provide. 



The final equation for the integrated DSPC, Eq. (16 1, 



reveals much about the eigenvalues of the percolating 
cluster. We can, for example, track down the origin 
of the 6 peaks in the density of states of the perco- 
lating cluster. Suppose for the moment that 6 c ,n(A) 
is continuous as a function of CI. Then the integrals 
over d6 Cj o(A) certainly do not generate any jumps in 
A perc (f2,c) due to the continuity in Cl. However, we 
know that a c * i a(A) = X) n a «^(A ~ z n) is discontinuous 
as a function of Cl since the z n depend on Cl, and in- 
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Figure 1: The integrated density of states at c = | and c = 2. 
The bottom figure shows an enlargement of the curve for c = 2 
around = 1. 



spcction of Eq. ( 16 ) shows that this leads to a jump in 
A perc (n) if the location z n of a jump moves from — oo 
to +oo when D, is increased infmitesimally (it does not 
lead to a jump if z n moves across or 1, as could be 
suspected at first sight, since the various contributions 
cancel in these cases). Eq. (13) shows that z n can and 
does indeed pass oo as is varied. While the peaks of 
the density of states of the finite clusters are located at 
those f2 for which z n — 0, the peaks for the percolating 
cluster are located where z n — oo. Since the zeros and 
the poles of z n necessarily alternate when Q is varied, it 
follows that the peaks in the percolating cluster lie dense 
if the peaks in the finite clusters lie dense. In this sense 
there is no mobility edge in the percolating cluster since 
isolated and thus localized eigenvalues exist throughout 
the whole range of < £1 < oo. 

Unfortunately, this argument only strictly holds if 



&c,n(A) is indeed continuous in f2. If it is not, cancel- 
lations might occur which could reduce (or even remove) 
the peaks from the percolating cluster. It is shown in [13] 
that indeed & Cj n(A) is not continuous but the argument 
presented there also shows that a complete removal of 
peaks would seem an extremely fortuitous cancellation. 
In order to check these results we have performed popu- 
lation dynamics simulations of the integrated DSPC ac- 
cording to Eq. (15). Fig. [I] shows the integrated DSPC 
directly at the critical point c = | and deep inside the 
percolating regime at c = 2. Some jumps are clearly ob- 
served for c = i, and they are located at the predicted 
positions. The most prominent ones occur at O = 1 



where the coefficient z = 

Q Q_ 

where z = 



oo, 



or at ft 



_ 3±VE 



n-i-i 



= oo. For c = 2 the integrated den- 
sity of states in Fig. [I] looks smooth. This is, however, 
not the case as the close-up in the bottom part of Fig. [T] 
reveals. The jump at fi = 1 which is very pronounced at 
c = | is also present at c = 2. 



Our work is based on ideas by Kurt Broderix who died 
on May 12, 2000. We thank Peter Miiller and Reimer 
Kiihn for valuable discussions. 
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